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CERTAIN  FINITE-DIFFERENCE  SCHEMES  FOR 
EQUATIONS  OF  THE  NONSTATIONARY  LAMINAR 
BOUNDARY  LAYER 

A.P.  Oskolkov 

Department  of  Applied  and  Computational 
Mathematics 

A number  of  stable  implicit  finite-difference  schemes  for 
the  solution  of  the  nonstationary  Navier-Stokes  w equations  is 
proposed  in  works  [1]  and  [2].  In  this  article  analogies  of  these 
schemes  for  equations  of  the  nonstationary  laminar  boundary  layer 
are  given;  it  is  shown  that  these  schemes  in  the  two-dimensional 
and  three-dimensional  case  are  decomposed  into  uniform  finite- 
difference  schemes;  the  unique  solvability  of  the  appearing  lin- 
ear algebraic  systems  of  equations  is  proven,  and  it  s4  is 
shown  that  these  systems  can  be  solved  by  means  of  a uniform 
trial  run. 


The  two-dimensional  nonstationary  flow  in  the  laminar  boun- 
dary layer  is  described  by  the  following  set  of  equations  and 
initial  and  boundary  conditions  [31,  [^]: 


du 


, J*u  dp  du  d9  „ 

fU ftf’— — — -t-  - — = "» 


du 


du 


dx 


(J  0^*  ~ dx  dx  0 y 

O^x^Xcoo  , Ossy-coo  f 0%t^T  , 

=0.'’!  =v(i,t),“(  ='('{». t) 

• 0 lyo  ' |x»0  '*  ' 

u (3c.y,t)  — u{x,t)  , 


(1) 

(2) 

(3) 

(4) 


1 


where  In  conformity  with  the  Bernoulli  equation 


Furthermore,  from  the  continuity  equation,  the  initial  con- 
dition for  u and  the  boundary  condition  for  v when  y = 0,  deter- 
mined Is  the  initial  condition  for  v 


The  solution  by  the  method  of  finite  differences  to  the  prob- 
lem formulated  by  dependences  (l)-(6)  is  virtually  impossible  due 
to  the  fact  that  the  variable  y is  changed  on  the  seml-lnf inlte 
interval  [0 , . In  work  [5]  the  author  goes  around  this  dif- 

ficulty in  that  by  means  of  the  Krokko  [Trans,  note:  spelling  not 
verified]  transform  f = 3:  , transforms  (l)-(6) 

into  the  initial-boundary  value  problem  for  function  W,  in  which 
the  variables  ? • 1 are  changed  in  the  finite  limits:  0^^«s;X  • 

. On  the  other  hand,  from  the  theory  of  the  boundary 
layer  it  is  n known  [31  that  u tends  to  u{cr,t^  when 

y — oo  exponentially,  and  therefore  the  is  boundary  condition  (4) 
can  be  considered  as  being  fulfilled  when  y = Y,  where  Y is  a 
sufficiently  large  finite  number,  whereupon  the  admissible  quan- 
tity Y can  be  defined  according  to  the  assigned  accuracy  with 
which  the  solution  is  sought  by  iterations.  Therefore,  subse- 
quently with  the  writing  out  of  the  finite-difference  scheme  for 
the  problem  (l)-(6)  we  will  assume  that  its  solution  is  sought 
in  the  domain  0<x<X  , and  . 

Let  us  divide  the  parallelepiped  * [O.Vj  * [0,T]  into 

elementary  cells  by  planes  x.“tAX,  where 

, and  we  assume 


r-ix 


Vy  S u ( tAX 


Then  by  analogy  with  the  Navler-Stokes  equations  [1],[2],  the 
system  of  equations  (1)  can  be  approximated  by  the  following  Im- 
plicit finite-difference  scheme: 


.f+i 


it 


At 


+ U 


u/:i-2uSu:'" 


l-  , / U7  — U.  U -2li  +u 

^ -iJ*i W'_l 


AX 


- '■t 

AX 


Ay  (Ay)* 

It  . 

= 0 , 


V 


r»^^4  , - 

= (7) 


(8) 


where 

1=1, 2, ...,L  , i=0,i.l N-1  . 

Added  to  equations  (7)  and  (8)  are  the  following  Initial  and 
boundary  conditions: 


II 

II 

% . t=0,{.. 

..,L 

(9) 

u^*‘=0 

to 

’ lO  t 1 

’ V ; ’ 

(10) 

il 

J=0.l,....N-i  . 

i-0,i 

.(II) 

Prom  equations  (7)  and  (8)  it  Is  clear  that  introduced  into 
them  are  values  of  function  U on  two  adjacent  verticals  and 

, and  values  of  function  v - only  on  one  vertical  s:.  , i=i,2,...,L  . 
By  using  three  of  the  boundary  conditions  (10),  we  obtain  the 
cyclic  process  for  the  find:!ng  of  the  unknowns 

each  ith  step  of  which  we  have  on  the  1th  vertical  at  any 
1=1,2,....,  N-l  a system  of  2(M-l)linear  algebraic  equations  with 
2(M-1)  unknowns,  and  these  equations  should  be  solved  by  taking 
Into  account  the  initial  and  boundary  conditions  (9)-(ll).  There- 
by the  problem  Is  reduced  to  the  subsequent  solution  of  the  uni- 
form finite-difference  problems  (7)-(ll)*  Let  us  show  that 
theorem  1 is  correct. 
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The  system  of  equations  (7)-(ll)  is  uniquely  solvable  at 
each  i={,2,...,L  and  each  i=0,{  . 

Since  the  equations  (7)-(ll)  arethe  system  of  linear  alge- 
braic equations  in  which  the  number  of  equations  is  equal  to  the 
number  of  unknowns,  then  for  proof  X of  the  theorem  it  is  enough 
to  show  ahX  that  the  corresponding  uniform  systems  have  only  a 
trivial  solution.  Further,  starting  from  the  uniform  initial 
condition  (9),  we  find  that  on  the  ith  vertical  when 

t = the  corresponding  (7),  (8)  uniform  system  has 

the  form 


at  (a^y 

i 

— ^ -^^=0. 

AX  <■}  aij 


(7a) 

(8a) 


where  f i = ,■  ■ ■ , H~l  • 

Added  to  these  equations  are  the  uniform  boundary  conditions 


, u^:‘=o , t=f.2....,L,  1=0,1 y-/, 

lO  ' tn 


(12) 


©•^*'=0  , 1-1,2,... ,L  , t = 0,f,.  .,N-l 


(13) 


It  is  known  [6]  that  the  system  (7a),  (12)  has  only  a 
trivial  solution  therefore  from  (8a)  and  (13)  it 

follows  that  = 0 = , J.=0,i,. ..  , -t  = 0,t,  . . . , N - i. 

Thereby  the  theorem  1 is  proven. 


Let  us  point  out  the  simple  method  of  solving  the  system 
(7)-(ll)-  For  this  let  us  rewrite  equations  (7),  (8)  in  the 
following  way: 
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u^'^  + C^  , 

=9'  , 

V A9  '■^  ay  »V 


where 


V (a^f  ^ ay  v)  ’ V ax  V (ay) 


-^-2  , f^=atF^)'+u^  + , g^.  = u^*‘ 

V (ay)*  V V ax  V ‘■''V 


Equations  (14)  together  with  the  boundary  conditions 


u^;'=o  , 

tO  > iM  <■ 


,;♦!  _ u/  '♦< 


the  initial  conditions  (9)  and  boundary  condition  can 

be  solved  on  each  1th  vertical,  i = <,2,-.  ,4  by  the  dispersion  method 
[7],  and  after  this  is  determined  from  (15)  with  the 

use  of  the  boundary  condition 


by  the  formula 


The  three-dimensional  nonstationary  flow  in  the  laminar 
boundary  layer,  circumfluent  a piece  of  plane  (x,  z),  is  des- 
cribed by  the  following  set  of  equations  and  initial  and  boun- 
dary conditions  [4]: 


flu  d u du  flu  , fl*u  flp 

dt  dx  dy  dz  fly  dx 

div  flW  flW  flw  d^W  flp 

_ +u  - -+tr--r— + w >)  — ’ 

dt  dx  dy  dz  oy  flz 

+ — — + = 0 , 0^x<X  , 0<i^Z  ,0*u«>o  ,0(tiTi 

dx  dy  fli  'i  ^ 


“Iv-o  ■ 


u{x,^,i,t)  — (J(x,2,t)  , W(x,y,z,t)— W(x,«,t),;y-«>  . (2i) 

Here  the  derivatives  of  pressure  ^-  = F(x,'i,t)  and 

Up  px  ooc 

are  considered  to  be  the  assigned  functions. 
Furthermore,  from  the  continuity  equation,  the  initial  condi- 
tions (18)  for  u and  W and  the  boundary  condition  for  v when 
y = 0,  it  is  possible,  u just  as  in  the  two-dimensional  case,  to 
define  the  initial  condition  for  v 


Having  in  mind  to  use  the  finite-difference  method  for  sol- 
ving the  problem  (17)-(22),  according  to  the  very  same  consider- 
ations as  in  the  two-dimensional  case,  we  replace  the  semi-infi- 
nite Interval  of  the  change  in  variable  y by  the  finite  inter- 
val [0,  Y],  where  Y is  a sufficiently  large  number  the  accessi- 
ble value  of  which  again  can  be  found  by  Iterations,  and,  there- 
by, we  examine  the  problem  which  describes  the  relations  (17)-(22) 
in  the  domain  0«x$X  , , and  we  replace  the 

boundary  conditions  (21),  respectively,  by  the  following: 


u -u(x,x,t)  , u/— , y— y . 

Let  us  divide  the  parallelpiped  5r=M’‘M’'[0.Z]x[0.T]  into 

the  elementary  cells  by  planes  'x.=iax  , • 

where  ax=^  , J=0./ M K»0.f R 

T 

~ > we  assume 

M 


1 


Similarly  determined  are  , lO"^^  and  so  on. 

By  analogy  with  the  Navler-Stokes  equations  [1]  and  [2],  the 
system  of  equations  (17)  can  be  approximated  by  the  following  im- 
plicit finite-difference  scheme: 


aI 


AX 


. _ j. 


'•■i*  AZ 


(^V) 


If,:. 


At  AX  r A^ 


At  ' V* 

AZ 


^ . (24) 

V f / 


O’.  . ” O’ 


where 


. ilf — ^ q.  , (25? 

AX  A A Z 


t-1,2,.  .,b  , J=l,2,...,M  , k=^/,?.  • • ,R  , l = 0,i,...,H-l  . 

Added  to  equations  (23)-(25)  are  the  following  initial  and 
boundary  condtlons: 


(26) 


(2T) 


(28) 


-®r 

V* 

..,L  , 1 = 0,1, 

, « 

=o,i.... 

0 , 

LCk 

» 10* 

LK 

4,  OK 

L^0,l 

,...,L,  A-l 

0,1,... 

, ^-o.t. 

t*l 

ttl  ■ .. . 

O^K 

- 

=j 

.,i  ,y- 

O^K 

t/" 

' 1^0 

1 

'A 

0,1,.. 

, 1^0, 
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(29) 


From  equations  (23)-(25)  It  is  clear  that  entering  Into  them 
are  values  of  functions  u,  w on  three  adjacent  verticals  with 
numbers  (t.^)  and  (i.n-i),  and  values  of  function  v only 

on  one  vertical  • Using  boundary  conditions  (28),  we  ob- 
tain the  cyclic  process  for  the  finding  of  i 

in  which  on  the  (<-.«) -th  vertical  at  any  i = X-i  we  have  a 

system  of  3(M-1)  linear  algebraic  equations  with  3(M-1)  unknowns, 
whereupon  these  equations  should  be  solved,  taking  the  initial  ■ 
and  boundary  conditions  (26)-(29)  into  account.  Thereby  the  prob- 
lem (17)-(22)  is  reduced  tothe  subsequent  solution  of  the  one- 
dimensional  finite-difference  problems  (23)-(29).  Just  as  in  the 
two-dimensional  case,  it  s is  shown  that  theorem  2 is  correct. 

The  system  of  equations  (23)-(29)  is  uniquely  solvable  at 
any  i=i,2,  ,L  , . 

The  system  (23)-(29)  permits  a simple  method  of  solution. 

For  proof  of  this,  let  us  rewrite  the  equations  (23)-(25)  in  the 
following  way: 


./  AMt  &L  t\  o I . at  . at  f 

yj'  ' V*  ^ /sx  ^ 42  (^' 

c«  ,f«  =_^Jf'''^u'  t— u'  u'  , 

f.^  =~AtG^*^+ u>^  .f u^  u>^*^  i- ~-uy-.^  uj-.^i‘ 

Ji^K  V*"  AX  v*^  ‘■■''^*  ~ V* 


Ai 


t Ay  t*i  ^y  „ 
q - ~~-u  , + — ? u>  . . 

?*■/*  4X  <--^,;»  4 2 Vi*'^ 


1(33) 


The  unknowns  determined  independ- 

ently from  each  other  from  equations  (30),  (31)  and  the  initial 
and  boundary  condtiions  (26)- (27)  by  means  of  the  one-dimen- 
sional dispersion , and  after  this  the  unknowns  are  deter- 

mined from  equation  (32)  and  the  second  of  the  boundary  con- 
ditions (27)  according  to  equation 


t 

‘■i*  vv 


-I- 

\ a: 


u: 


iH 


az 


By  applying  the  h known  sufficient  conditions  of  correct- 
ness of  the  one-dimensional  dispersion  [7],  and  using  the  fact 
that  according  to  the  physical  meaning  of  the  problem 

j it  is  easy  to  sia  establish  that  the  solutions 
to  equations  (30)-(33)  are  a stable  with  respect  to  the  kh  com- 
putational errors  if  the  following  condition  is  fulfilled; 


+ 


min  15-/ >0  . 

ijKt  ‘Z* 


In  the  two-dimensional  case  the  condition  of  correctness  of 
the  dispersion  for  the  system  of  equations  (14)-(16)  has  the  form 

— + min.  . 

Ay  ijl 


I 
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